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ABSTRACT 

An attempt was made to forecast the occurrence of llontercy procipita= 
tion using only 500=mb data. The use of tke wind spoed and direction and 
the 24-hour wind direction change at weather ships Papa and ilovember gave 
no significant results for the poriod of December 1955 to March 1956, 

Using West Coast stations, a probability forecasting graph was 
plotted from the height differonce between Santa Liaria and Oakland and 
the 500-mb temperature at Oakland. The results on an independent sample 

sof 151 casos resulted in 75% correct forecasts with a skill score of 0.29. 

With a sample of 72 tho 2-parameter forecast was 68, corroct with a 
skill score of -0, 08, When a third parameter, 500-mb wind direction at 
Oakland, was added, tho accuracy increasod to 73% corroct with a skill 
score of 0.11. For this set of data neither method gave significant 
results. 

Although these results do not show the accuracy that other objective 
methods generally havo, these rosults show the possibility of forecasting 
surface precipitation from 500-mb data, 

Iho study of the winter time precipitation at Monterey was undortaken 
in the spring of 1957 at the U. S. Naval Postgraduate School with the 
purpose of forecasting the occurrence of lontorey precipitation fram 
o00emb data. The author wishes to express his appreciation to Professor 
We De Duthie for his valuable guidance and heartening enoouragemont. The 
author is also grateful to the Engineering Department of Dol Monte 


Properties Company for freo aooess to their weather obsorvational data, 
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De Introduction 

Objective forecasts of two t:pes are presently used to produce 
fastor, moro accurate forecasts and to supplement the technique of fore= 
casters who are less experionced in the idios:ncracics of spocific 
localitics. The user of Vernon's method [1] for prodictinzs precipitation 
at San Francisco orjectively types weather situations and then uses sur- 
face information to find tho forecast of "rain" or “no rain." 271 in 
his forecast of precipitation for Washington, D. 2. [2] and the TVA Basin 

[3] , developed a method for forcoasting the probability of precipitation 
using as many parameters as are necessary and pertinent. On the West 
Coast, Jorgensen [4] applied this method to forecasting the occurrence of 
rain in the San Joaquin Valley of central California, and Thompson [5] 
similarily quantitatively forecast precipitation for Los Angeles. 

Ihe probability forecast has a definite advantare over the 
categorical forecast. То be useful the categorical forecast must be 
botter than, or nearly as good as the subjective forecast. + Although an 
experionced Torecaster may have a foclins of how accurato his subjective 
forecast is, few if any can or have ever tried to express their confidonce 
in the form of a percentage. On the other hand even a probability fore= 
casting method which does not quite measure up to subjoctivo forecasts in 
the categorical field may accurately indicate the probability of rain in 
every instance. The categorical forocast may be thought of as having the 
50/0 frequency line as the divider between "rain" and "no rain." This 
50% line may not always be the deciding line for attempting an operatione 
One must decide in each case by the cost of the operation in money or 


human lifo, the gain to be expected, and the loss to be expected from 





the failure of the weather clement. The maximum or minimum probability of 
rain for the operation under consideration may be calculated. This process 
of acting on calculated risks, though universally used for a long time in 
other economic practices, had to wait for probability forecasting to be 
applied to the weather. The U. 5. Weather Bureau has found their experi- 
ment in issuing probability forecasts successful. Even the Sunday afternoon 
picnicker wants to know how sure it is to be a clear day. 

Since the Numberical Prediction Center in Washington, De Ce produces 
a much better 500 than 1000-mb prognostic chart, the present development 
of a probability forecasting technique for the occurrence of rain in 
Lionterey used only 500=mb datas The method itself was developed using 
actual radiosonde data, but it was expected that the forecast might be 
extended using the numerical prognosis at 500 mbs. This present study ів 
unique in that all other methods in the literature depond to a certain 
extent on surface data. Even Jorgensen [4] used surfaco data to type 
weather situations before forecasting with 700-mb parameters. ihe аГогс- 
mentioned methods relied solely on the 700-mb level for their upper-air 
data. low that there is a sufficient history of accurate 500-mb data, 
relianco on this lower level for all uppereair data is not necessary. 

The use of the raw data from radiosonde reports has the advantage of 
being more objective than the subjectively analyzed data. The use of 
teletypo data also allows the forecaster to solve his problem a couple of 
hours sooner 17171 if he needed an analyzed chart. The disadvantages of 
using this raw data are the possibility of unrepresentative data and 
missing or garbled reports. In spite of the fact that the data at 500 mbs 


is much more representative of the area than that below, and that this 





data is recorded and transmitted much moro reliably than beforo, forecasts 
on a “perfect 500=mb prognosis” probably would be as sood as those using 
actual datae 

To verify Vernon's hypothesis [1] that his method could bo used 
satisfactorily without modification within a 150 mile radius of San 
Francisco, llacLeod ka and Jepson [8] testod this mothod at Шопбогоу, 
Although the results azsreed quite well with Vernon's results [1] و‎ 6 
should not assume that the pattern for iiontorey precipitation is greatly 
similar to that of San Francisco. Even on tho Lontoroy Poninsula there 
may be considerable variation. The precipitation measurements for this 
research wero taken at Pebble Beach. In gencral ionterey precipitation 
occurs at the samo time in somewhat larger amounts. Tho fact that this 
method was specifically developed for Pebble Beach precipitation should 
not hinder its use anywhero on the peninsula, ospocially in view of the 


studies of MacLeod [7] and Jepson [8] . 





2. Preliminary Studios on the 1955--56 Rainy Season 

In this initial work precipitation measurements from the laval Air 
Facility for December 1955 through March 1956 were used. In this period 
of 122 days 46 rainy days were recorded. 

In order to test the local forecasting rule that a SW wind in lionterey 
is followed by rain the next day, the concurrence of those two phenomena 
was tabulated. Although precipitation and SW winds individually occurred 
46 times, they coincided in the proper manner on only 14 occasions or 30% 
of tho times During this particular year the rule had no significance, for 
rain itself averaged a 39% chance of oocurrence. This short sidetracr 
was the only deviation from the use of 500-mb data as prodictorse 

Attempts wore mado to correlato tho 500-mb radiosonde data from 
weather ships "Papa" and "November" with rain at the Naval Air Facility. 
The first test involved the correlation of a 24-hour back or veer of the 
wind and subsequent precipitation in Monterey. Table 1 shows the results 


of this experiment. The total number-of cases for each column is less than 


No of days later Papa „.ovember 

rain occurred Back Veer Back Veor 
0 16 55% 24111 13 39% 8 42% 
1 9 له‎ 4 46 16 48 5 6 
2 12 41 7 54 15 45 77 “94 
2 15 45 7 54 15 39 5 26 
4 11 38 85 72 14 42 6 32 
5 11 98 8 62 17 51 4 21 

Total No. Cases 29 13 33 19 

Table 1 


the sum of each column because each variation in wind direction was often 
the forerunner of several rainy days. Montorey precipitation generally 
persists for a few dayse The verification percentages are, for each type, 


the specified number of days after the indicated change of wind direction. 
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Even from this small sample it can be seen that this chanze docs not in 
itself provide a significant forecasting tool. Actually the best rosults 
were found for four to five days after veering happened. Even though 
weather ship Papa is a considerable ‘distance upstream from l.ontéerey, it 
is hardly conceivable that events happening four or five days previously 
could be more important than more recent events. Actually on the basis of 
the size of the sample there is little significant difference between eash 
day. There does seem to be some significance in the fact that the verifi- 
cation percentazes for Papa increase for veering winds and decrease for 
backing winds as the length of the forecast period increasose Ihe 
observations for weather station November seem to be more random with less 
differentiation betweon veering and backing. 

An attempt was made to show the effects of wind direction and speed at 
Papa and lovembor and subsequent rain in Monterey. Graphs for the occur- 
rence and non-occurrence of rain after zero, one and two days were plotted 
for November and Papas -igures 1 and 2 show these plots. Tho length of 
the forecast period showed Ne Significance. Using the mean valucs for these 
days, Figuro 3 was plotted to show the difference between days of rain and 
no rain. Hore again the two plots did not differ significantly. Although 
the frequencies differed, the maxima and minima for the rain and no-rain 
casos came at the same wind speeds or directionse 

From these preliminary studies it was concluded that although reo- 
ports from Papa and liovembor might be connected with subsequent rain in 
llonvercy, more usable information might be obtained using stations closer 


to Monterey. 














= bh — 

— | - 
8 Г 

т”? 17 | 
— - | 1 ЕЙ 
| چم‎ | 
اسم — وھ‎ > 
AOS 


Е lot 


i]!‏ = ,سيم 
mas‏ | 


im 4 























Je Use of West Coast Reports 

Since the failure to find any corrolation between the 500-mb reports 
of the weather ships and the occurrence of Lionteroy precipitation could 
be blamed to a large measure on the distance betweon these stations and 
llonterey, the next effort was devoted to the oorrolation of West Coast 


o00-mb reports апа Monterey rain. Table 2 lists the twelve winter months 


1949 January February March December 
1950 January Fobruary March 
1951 à; March December 
1952 January February Maroh 

Table 2 


which were used for this study. All upper-air data were taken from the 
Synoptic Series [6] e The precipitation reports taken by the Del Monte 
Properties Company at Pebble Beach were used because they were the best 
available source of data. The available California Summaries were 
incomplete, Rain measurements from the Naval Air Facility differed 
slightly from those taken at Pebble Beach; therefore all the results of 
this study should only apply to Pebble Beach. 

In the Pebble Beach data there were some days in which no record was 
entered. These days, along with the days in which the pertinent informa- 
tion was missing, were not used. In spite of these omissions, approxi- 
mately 350 cases were used for cach study. Brier [2] used 271 cases for 
precipitation studies in the ¿=month period from December through February. 
Thompson [5] employed approximately 360 cases for his precipitation study 
covering the six months from October through Larch. On the basis of 
these past studies, the present sample seemed adequate. 

Three stations were chosen so that the lateral variation of the 
parameters could be measured. Since nearly all storm activity comes 


9 





from northwost of Monterey, Santa AT and Oakland were chosen to 
bracket Monterey, and Medford was selected to indicate the approach of 
weather-producing phenomena from tho north. "hen à single station was 
used, the Oakland reports were naturally used because this station is 
both near and north of liontereye 

The 500=mb parameters recorded for each station were the height of 
this surface, the temperature and the wind direction and speede Prelimi- 
nary scatter diagrams using A Horry САК у5 А Нодк-мер 9nd the plot of 
A TK OAK against TT oak were used to decide the most favorable length 
for the forecaste Beside each point the number of days until precipitation 
occurred was plotted. When precipitation persisted for several days, the 
number of days till the termination of precipitation was also plotted 
since the parameters might be construed as forecasting rain for any of 
these days. Figure A is an example of one of these plots. Irom these 
scatter diagrams the best forecast was determined to bo for the same or 
the following day. The radiosonde reports published in the Synoptic 
Series [6] are recorded for 03002, which is 1900 PST of tho previous day. 
The Pebble Beech rain gauge was road daily at 0815 PST, and therefore the 
length of time until the beginning of the forecast period was 13.- hours. 
At the time of receipt of teletype data this forecast would be for 
approximately ten hours. Thompson [5] estimated a six-hour forecast 
after map time. Brier [1] forecast for a 24-hour period starting six 
hours after a local surface observation. Thus if significant results 
can be obtained in this study a considerably longer forecast period may 
be used, even in cases in which missing data requires the partial analysis 


of the 500-mb chart. In order to overreach the 24-hour prognosis of 
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Vernon [1] the use of a good prognostio 500=mb chart would be necessary. 

In the sraphical integration technique suggested by Brior [2] , any 
number of variables may be combined into one forecast. Two parameters are 
plotted against ono another, and each point is markod appropriately for tha 
occurrence or non=ocourrence of the meteorolorical evont. This scatter 
diagram is then analyzed into isopleths of percentage occurrences The 
values of these isopleths aro used as a new variable derived from the 
original pair in the construction of a new scatter diagrame, Іп this step» 
wise manner any numbor of variables. may be combined into one forecast. At 
each stage the derived variable should represent a more definitive forecast. 
The procedure seems very simple, but the accuracy of the forecast depends 
on the selection of thé most significant variables and the proper analysis 
of the scatter diagrams. | | 

In the present application of Brier's technique [5] variables were 
selected on the basis of thoir availability as well as their forocast 
significance, so that a reasonably good forecast could be made Prom tele- 
type data with partially missing reports. Dew=-point temperatures wore 
omitted because they wero so often missinge Since wind direction and speed 
were frequently absent, they were omittod from the "first approximation" 
to tho forecast. Only the heights and the temperatures at the 500-mb 


level were used initially. 
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Le гр Corrolation of liontorey Procipitation with Bs سوم‎ and 
OAK MER 
Figure 5 is the scattor diagram of C He rx OAK plotted against 
Glo k.Mpg* Crosses indicate no subsequent procipitation at Pobble Beach, 
and numbered ciroles indicate tho amount of rainfall. This diagram was 
then divided into aroas having ten or more obsorvationse The fractional 


2 test 


and percontase occurrences of rain were indicated in each aroa. AX 
[9] was run to show that the distribution of points was significantly 


different from chance. The formula for evaluating x is 


x2 = $ lo el? 


where "о" is the observed frequency and "e" is tho chance frequency. The 
level of significance of this sum is then found from a table of X? values 
tabulated against the numbor of degrees of freedom. Because the X? test 
is only meaningful for valuos of e greater than five, adjacent blocks 
having similar composition wore combined as illustrated in Table 3. The 
rosults of this test show that the distribution of points within the cells 
differs significantly from chance well beyond the 0.001 confidence levele 
This agreed with the results of Brier [2] a 

The X? test is a necessary but not sufficient condition for the oon- 
struction of a proper percentace forecasting chart. In this test no 
account can be taken of the juxtaposition of blocks of widely differing 
percontages of occurroncee In this instance the isolines could not be 


plausibly dram, and thereforo the diagram had to be abandoned, 
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Block Nose Noe of Cases Obsorved Expected (se) 





ec ы 29 3 10.4 5.26 
4 & 5 17 9 6.1 1.38 
6 6 7 20 14 7.2 6.42 
8, 9,810 37 0 13.3 13.30 
Me 12 22 7 7.9 0.12 
13,14 & 20 32 3 11.5 6.28 
15 & 16 20 9 7-2 0.45 
17,18 & 19 35 15 12.6 0.46 
22 & 23 27 12 9.7 0.55 
25 & 26. 26 , 18 9.3 8.14 
24,25 & 28 33 17 11.8 2.29 

24.66 


lO degrees of freedom X^ = 29.6 at the 0.001 level 


Table 3 
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5. Correlation of lionterey Preoipitation with А Н and TTo 


5МХ-ОАК АК 

From a marginal analysis of Figure 5 it was determinod that the ocour- 
rence of Monterey precipitation was more highly correlated with A How OAK 
than with Alo Therefore the former variable was retained and 
plotted against TT oaks This plot had the advantage of requiring data orly 
from two stations instead of threc. In the division of the scatter diagram 
into blocks an attempt was made to make the blocks as nearly square as 
possiblo. In this manner the peroentage of occurrence for each block can 
be said to represont that of the center with the least error. Figures 6 
and 7 show the scattor diagram and the block diagram, respectively. 
Although the X% test of the data, as shom in Table 4, did not show signifi- 
oance quite to the 0.001 level, analysis into isopleths of percentage as 


illustrated in Figure 8 was easily accomplished. 


Block Nos. Ho. of Cases Observed Ехрос+ей (о-е)2/ 

1'& 2 23 13 8.3 2.66 

3 7&8 36 18 13.0 1.92 

4 15 11 5.4 5.69 

9 & 12 22 11 7.9 1428 

5, 15614 35 19 12.6 3.26 

10 & 11 27 9 9.7 0.05 
6 & 15 28 5 10.1 2.58 

16 17 4 6.1 0.72 
12 & 17 21 7 7.6 0.04 
20, 22 & 23 27 10 9.7 „01 
18 15 0 5,4 5.41 
19: 18 3 6.4 „81 
21, 24 4 25 33 4 11.9 5.24 
50.60 


12 degrees of freodom X? = 32,4 at the 0.001 level 
Table 4 
This prediction graph was tested against data for January, Fobruary 
and larch of 1953. Since the data from Pebble Beach were not recorded on 


Sundays in this year, no forecasts for Sunday and Monday were made when 
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the reading for llonday showed a two=day accumulation. Since this sample 

ОГ 79 was too small to test the breakdown of each ten porconta;o line, 

this data was only tested on a categorical forecast using 50% and groater 

as the area of the graph for forocasting rain. Table 5 shows that the fore- 


cast verified 815 with & skill soore of 0,45, 


Forecast 
Rain No Rain Total 
Rain 8 7 15 
Observed No Rain T 97 64 
Total 15 64 79 
Table 5 


The skill score is defined as: 


0 - Ez 
I - Eg 
where C & the number of correct forecasts, 


So 2 





E¿* the number of forecasts expected to be correct due to chance, 

I = the total number of forecasts. 
So has a value of unity when all forecasts are correct, and a value of 
zero when the number of correct forecasts is equal to the number expected 
by chance. Negative skill scores arc possible. Eg 15 definod as follows: 

Ес = Fp(R/T) Fp (R/T) 

where Fy = the number of rain forecasts, 

Fars the number of no rain forecasts, 

R = the number of observed rain periods, and 

ЦВ = the number of observed no rain poriods. 

In order to test the forecasts against a skilled forecaster who has 


all the current data, these forecasts were compared with those the U. Se 
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Weathor Purcau published in the lontorey Peninsula ¡¡crald. Since this 
paper does not publish on Sundays, no forecasts wero available for llondayse 


From Table 6 tho verification was 865 with a skill score of 0.67. The 


skill score was significantly better for this subjective forecast than for 


the first approximation of the objective forecast. 


Forecast 
Rain No Rain Total 
Rain 15 6 ЕП 
Observed No Rain 4 48 92 
Total 19 54 75 
Tablo 6 


Data from December 1946, liarch 1947 and December 1948 were addod to 


tho first tost data and vorification for the six-month sample is shown in 


4 


Tablo 8. The forecasts wero 75% correot with a skill scoro of 0.29. The 


last three months tested were not compared with U., Se Woather Bureau 


forccasts, 
Forocast 
Rain Хо Rain Total 
Rain 15 18 33 
Observed Мо Rain 19 99 118 
Total 34 117 151 
Table 7 


Although these results were somewhat discouraging, they actually are 
bettcr than those obtained by using the original data. from Table 9 


there was 65% verification with a skill score of 0.22. 
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Forecast 


Rain llo Rain Total 
Rain 54 67 121 
Observed No Rain 50 161 211 
Total 104 228 332 
Table 8 


Table 9 shows the breakdown of the original data and the independent 
data for each ten=-percent interval. The class mark, plotted as the ex- 


pected value, is the mid-point of the interval. 


Original Data Independent Data 
Expected No. of Cases % Rain Но. of Casos ^ Rain 
457 58 10 17 18 
145 43 12 7 29 
245 51 21 7 0 
ae 48 21 12 25 
42 58 45 10 50 
52 40 5 9 20 
648 57 70 A, 20 
745 25 52 6 50 
Table 9 


The failure to verify each area even for the original data must be 
due to the procedure of constructing the isopleths of percentage. The 
fact that the block diagram, Figure 7, was made from the scattor diagram, 
Figure 6, by constructing squares containing approximately the same number 
of observations automatically gives a larger weighting factor to the 
larger squares which have the fewer observations. Possibly the lines 
could be drawn more accurately on a chart which had the scales distorted 


so that each unit of area would have about the same number of observationse 
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6. Preparation of a 3=Paraneter Torecast 

In order to introduce new parameters not involving temperatures and 
heights, the first approximation to the percentage ocourronce of rain was 
plotted separately against the wind speed and the wind direction at 
Oakland. The graph of wind speed, ¡"“igurc 9, was abandoned with a Six- 
month plot because the wind speod did not seem to be producing a further 
separation of "rain" and "no rain." The 500=mb wind direction at Oakland 
soemod to give moro significant results. This scattor diarram, Figure 10, 
was analyzed into 15 blocks having 14 or more observations eache 
rigure ll shows this block diagram, which was converted into the forecast- 
ing graph of liguro 12, 

Of the independent data previously used, only the chronologically 
earlier data iricluded 500=mb windse This three month period forecast 
correctly 68% with a skill score of -0.08 using the 2-paramctor forecast 
and 73% correot with a skill score of 0.11 for the 3eparameter forecast. 
The second prediction graph improved the initial forecast only slightly. 
With this small sample and these low skill scores one can only conclude 
that if it improves the forecast, it improves it only slightly. Table 10 


shows the breakdown of the independent data for the J=parameter forecast. 


Fores ast 
Rain No Rain Total 
Rain 4 19 25 
Observed No Rain 8 41 49 
Potal 12 60 72 
Table 10 


It is interesting to note that the arca in Figure ll right next to 
the area of 8475 occurrence of rain did not have a single instance of rain 
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in the eight observations both in the dependent and independent data. 
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